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ABSTRACT 



We present deep (V - 28.0) BV photometry obtained with the wide field channel of the Advanced Camera for Surveys on board HST 
for four M3 1 globular clusters that were identified as candidate intermediate-age (age~ 1-9 Gyr) by various authors, based on their 
integrated spectra and/or broad/intermediate-band colors. Two of them (B292 and B350) display an obvious blue horizontal branch, 
indicating that they are as old as the oldest Galactic globulars. On the other hand, for the other two (B058 and B337), which display 
red horizontal branches, it was not possible either to confirm or disconfirm the age estimate from integrated spectra. The analysis of 
the distribution in the spectral indices Mg2 and H/3 of the M31 and Milky Way clusters whose horizontal branch can be classified 
as red or blue based on existing CMDs, strongly suggests that classical age diagnostics from integrated spectra may be significantly 
influenced by the HB morphology of the clusters and can lead to erroneous age-classifications. We also provide the CMD for another 
two clusters that fall into the field of the main targets, B336, an old and metal-poor globular with a significant population of RR Lyrae 
variables, and the newly discovered B531, a cluster with a very red red giant branch. 
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1. Introduction 

Nearly all types of galaxies contain globular clusters (GC), from 
dwarfs to giants and from the earliest to the latest types. Our 
most in-depth understanding of individual GCs comes from 
studies of the Milky Way (MW; see, for example, Dotter et 
al.|2010 and references therein). 

Among the external galaxies, M31 plays a special role in 
our studying of GCs. It is our nearest bright spiral galaxy neigh- 
bor, it is the most prominent member in the Local Group (LG), 
and it hosts the largest population of GCs: 544 confirmed GCs 
are listed in the latest version of the Revised Bologna Catalog] 
(RBC v4.0, Galleti et al. USS . At th e distance of M31 (~ 783 
kpc, from McConnachie et al. 12005b . GCs have almost stellar 
appearance (10 pc correspond to ~ 2.6"), hence the diagnostics 
based on integrated light that are used for more distant galaxies 
can be homogeneously applied here. On the other hand, M31 is 
also close enough that individual stars in GCs can be resolved 
and measured with the Hubble Space Telescope (HST). In this 



Send offprint requests to: sibilla.perina2@unibo.it 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program GO-10631 [P.I.: T. 
Puzia]. 

** Photometric catalogs are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-t-A/ and at 
http://www.bo.astro.it/M31/hstcatalog/ 
r http : //www . bo . astro . it/M3 1/ 



framework, the GCs system of M3 1 is a fundamental test bed for 
checking the consistency between the cluster parameters (in par- 
ticular age and metallicity) derived from the integrated light and 
those obtained from the analysis of the color-magnitude diagram 
(CMD) of individual cluster stars. 

The most powerful method to determine cluster ages relies 
on the accurate location of the main-sequence turn-off (MSTO) 
point (Renzini & Fusi Pecci 1988). Unfortunately, this method 
is limited to the nearest GCs where individual stars can be re- 
solved and measured down to at least ^1-2 magnitudes fainter 
than the MSTO. At the distance of M31 the task is made very 
challenging owing to the effects of crowding and the intrinsic 
faintness of the feature. At present only one of M31's GC has a 
CMD down to the TO, which has been achieved only with ~3.5 
days of HST/ACS integration time (B379; Brown et al.|2004f] 

In general, typical CMDs of M3 1 GCs obtained by HST data 
barely reach the base of the red giant branch (RGB), includ- 
ing most horizontal branch (HB) stars, except the bluest/faintest 
ones, when present (Ajhar et al. 1996; Rich et al. 1996, Fusi 
Pecci et al. [T9961 H olland et al. 1997; Jablonka et al . 120001 
Meylan et al. 200 H Rich et al. 120051 Per ina et al. I2009bt 
Huxoret al. 120041 [20051 120081 Galleti et al. 120061 Mackey et 
al. 120061 12007I I. In this magnitude range the HB morphology 
can be used as a rough age indicator, because the presence of 
RR Lyrae variables and/or blue HB stars implies low-mass pro- 
genitors, hence old ages. In a recent analysis, Dotter et al.( 2010) 

2 It is interesting to note that B379 was found by Brown et al. d2004l > 
to be 2-3 Gyr younger than Galactic GCs of the same metallicity, with 
an absolute age of lO^ 25 Gyr. Re-analyzing the same data and using 
different stellar models, Ma et al. d2010l l estimated an age of 11.0 ± 
1.5 Gyr. 
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conclude that all Galactic GCs with a blue horizontal branch 
have ages ^ 12.0 Gyr, on a scale in which the oldest GC is 
- 13.5 gyr old. 

To date, there are relatively few M31 GCs with a CMD (44). 
As a consequence, our overall characterization of the GC system 
of that galaxy mainly relies on the analysis of integrated colors 
an d spectra (se e Galleti et al. 120041 Puzia et al. 12005 1 Caldwell et 
al. 120091 1201 II and references therein). A widely used technique 
is based on Lick spectral indices (Burstein et al. 1984; Faber et 
al. 119851 Worthey[l994). Typically, indices that are mainly sensi- 
tive to age (Hp, H y ) are used in combination with others that are 
more sensitive to metallicity (Mgb, [MgFe]), attempting to break 
the well-known age-metallicity degeneracy. At fixed metallicity, 
stronger Balmer lines indicate a hotter MSTO, hence younger 
ages (see for example Puzia et al. 2005 and Caldwell et al. 2009 
1201 11 ). Another widely used method to age-date extragalactic 
star clusters relies on multi-band photometry: the age can be 
derived from the spectral energy distributions (SED) measured 
in broad-band photometric systems, by comparing them with 
simple stellar population (SSFp) synthesis models (e.g. Jiang et 
al. 120031 Fan e t al. E 556] 120081 l20T0l Ma et al. 125571 125091 
and Wang et al. 1201 Oi l. Clearly, photometry -based technique are 
prone to errors owing to the uncertainty in the interstellar extinc- 
tion, while spectral indices are virtually free from this effect. 

In a recent analysis adopting the most up-to-date theoretical 
tools and methodology, Puzia et al. (2005; P05 hereafter) derived 
spectroscopic ages, metallicities and [a/Fe] ratios for 70 GCs in 
M31 based on Lick indices. Within this sample, the authors find 
a population of ~20 GCs with ages between ~ 5 and 9 Gyr and 
mean metallicity of [Z/H] ~ -0.6. Independently, Burstein et 
al. ((2004) and Beasley et al. (2005) also claimed to have found 
two (B232 and B311) and six (B126, B292, B30L B337, NB16, 
NB67) candidate intermediate-age GCs in M310 respectively, 
in common with Puzia's sample. The presence of a population 
of intermediate-age GCs in M31 was also supported by other 
spectrophotometric analyses (Fan et al. 120101 Wang et al. 120101 
and reference therein; see Table EJ. 

On the other hand, Strader et al. (2009) showed that three 
of the candidate intermediate-age M31 GCs studied by P05, 
Burstein et al. (120041 ), and Beasley et al. ( 120051 ) have M/L v and 
V - K colors typical of old GCs. Twelve of the spectroscopically 
identified intermediate-age GC candidates were also detected by 
Rey et al. (2007 ) in far-ultraviolet (FUV) and/or near-ultraviolet 
(NUV) GALEX images. From the comparison of their UV pho- 
tometry with theoretical models these authors concluded that 
most of these spectroscopically identified intermediate-age clus- 
ters may in fact be classical old GCs with many blue HB stars. 
These hot stars would contribute to enhance the UV flux as well 
as the strength of Balmer lines, thus mimicking younger ages. 
A possible example is provided by the cluster B311, reported to 
be 5 Gyr old by Burstein et al. (2004), that reveals an extended 
HB in the CMD presented by Rich et al. (120051) . clearly not con- 
sistent with such a young age. Therefore, deep HST imaging 
of intermediate-age GCs candidates is crucial to check if they 
are genuinely intermediate-age or, instead, old clusters with blue 
HB. The presence of bright intermediate-age clusters in M31, if 

3 A simple stellar population is an ideal population of stars that have 
the same age and chemical composition (Renzini & Fusi Pecci 1988 ). 

4 It is interesting to note that Beasley et al. ([2005 1 find remarkably 
different age estimates from the same data, depending on the set of 
theoretical models they compare with. In particular, ages obtained by 
comparison with the Bruzual & Chariot (2003) models are systemat- 
ically lower than those derived from the comparison with Thomas et 
al. J2003t ones (see Table|2] below). 




Fig. 1. Location of the 4 target fields (open squares with the size 
of an ACS-WFC field) projected against the body of M31. Note 
the different orientations of the F606W and F435W exposures. 



confirmed, would be especially interesting because these objects 
are not observed in our own Galaxy, possibly implying a signifi- 
cant difference in the star- (or, at least, cluster) formation history 
between the two galaxies. 

In HST cycle 14 deep ACS/WFC imaging was acquired 
of four M31 GCs, which were identified by P05 as candidate 
intermediate-age cluster. In this paper we present the CMDs of 
the four main targets of this observational program, discussing 
the compatibility with the parameters obtained by P05 from the 
integrated spectroscopy (see Tables [2] and [3]). We also present 
the CMDs for two other GCs that were in the target fields in both 
passbands: B336 and the newly detected cluster B531. Finally, 
we provide firmer classification for three additional objects listed 
in the RBC that were imaged only in one passband. 

In Sect. [2] we describe the HST/ACS data, the adopted reduc- 
tion procedure, the photometric uncertainties and the complete- 
ness of the data. Sect. [3] is devoted to describe the CMDs of the 
individual clusters, the field-decontamination procedure, and the 
method we used to estimate the metallicity, reddening, and dis- 
tance. In Sect. |4]we describe the method for searching variable 
stars and the results for the individual clusters. Finally, Sect.[5]is 
dedicated to the discussion of the results and conclusions. 
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survey. 



ID 



RA(J2000) Dec(J2000) V„„ E(B-V)* 



B058-G119 
B292-G010 
B337-G068 
B350-G162 



00 41 53.01 
00 36 16.59 
00 40 48.45 
00 42 28.33 



+40 47 08.6 
+40 58 26.6 
+42 12 11.5 
+40 24 50.2 



14.97 
16.99 
16.73 
16.74 



0.13+0.01" 

0.13* 
0.06+0.02 
0.10+0.02 



The integrated apparent magnitudes of the clusters (V;,„) are taken 
from RBC V 4.0 (Galleti et al l2004l . * from Fan et al. f2008T >. ' Rich 
et al. J2005t find E(B-V)=0. 12 for this cluster. * from Caldwell et 
al. J201 li . These authors provided reddening estimates also for the 
other clusters, in particular E(B-V)= 0.15, 0.10, and 0.12, for B058, 
B337, and B350, respectively (uncertainties not reported). 
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Fig. 2. F606W images of the four candidate intermediate-age 
clusters that were the main target of the observations. The size of 
each image is 10"x 10". 10" corresponds to 38 pc at the distance 
offvDl (D = 783 kpc). 



2. Observations and data reduction 

The four main targets are the clusters B058, B292, B337, and 
B350. They were selected from Table A. 3 of P05 as candidate 
intermediate-age GCs.The observations were performed under 
program GO-10631 (P.I.: Thomas Puzia), using the ACS on 
board HST, and the filters F435 W (similar to Johnson B; four im- 
ages with a total exposition time of 4800 sec for B058, and 4900 
sec for the other three targets) and F606W (similar to Johnson 
V; six images with a total exposition time of 7500 sec for each 
target). Fig. [T] shows the location of the four target fields on the 
body of M31. In Table [T] we give the basic parameters of the 
target clusters. Their images in the F606W band are shown in 
Fig.0 
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Fig. 3. Stamp-images of other objects listed in the RBC that were 
included in the analyzed WFC images (B049, B336, SK087B, 
SK088B, SK089B) plus the newly discovered cluster B531. The 
size of each image is 10"x 10". 



The pointings were arranged to contain as many clusters or 
candidate clusters as possible because other confirmed or can- 
didate GCs are present in the proximity of the targets. This op- 
timized the scientific output of the images. For the fields with 
multiple objects, orientation requirements were set to guaran- 
tee the best coverage of all targets, and additional orientational 
constraints were added to the B292 pointing to avoid a bright 
star in the field of view. The field of B058 also contains B049 
(already studied in Perina et al. 2009b); the field of B337 also 
contains B336 and in the field of B292 we discovered a new 
cluster that was never detected before (dubbed B531, according 
to the RBC nomenclature). For these clusters it was possible to 
obtain meaningful CMD that are presented and discussed be- 
low. In the field of B058 we identified another three candidate 
GCs listed in the RBC (SK087B, SK088B, SK089B). These are 
faint and very compact objects for which we provide only a re- 
classification based on visual inspection of the ACS images. In 
Tableware listed some basic parameters of these additional ob- 
jects, and their images in the F606W band are shown in Fig. [3] 

The data reduction was performed in the same way as in 
Perina et al. (2009b ; P09b hereafter) on the pre-reduced im- 
ages prov ided by STScI, using the ACS module of DOLPHOT0 
(Dolphin 2000a), a point-spread function fitting package specif- 
ically devoted to stellar photometry from HST images. The 
package identifies the sources above a fixed flux threshold on 
a stacked image and performs the photometry on individual 
frames, accounts for the hot-pixel and cosmic-ray masking in- 
formation attached to the observational material, automatically 
applies the correction for the charge transfer efficiency (CTE, 
Dolphin 2000b ) and transforms instrumental magnitude to the 
VEGAMAG and standard B VI system using the transformations 
by Sirianni et al. d20051 >. 

We fixed the threshold for the search of sources on the im- 
ages at 3<t above the background. DOLPHOT provides as output 
the magnitudes and positions of the detected sources as well as 
a number of quality parameters for a suitable sample selection, 
with regard to of the actual scientific objective one has in mind. 

Owing to the adopted peculiar pointing strategy described 
above, we were forced to reduce the F435W and F606W frames 
of each cluster independently. The output positions were cor- 



See http ://purcell . as . arizona. edu/dolphot/ 
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Table 2. Ages from integrated spectroscopy/photometry for the main targets. 



ID 






AGE [Gyr]" 








P05 


Be05+BC03 


Be05+T03 


W10 


F10 


B058-G119 
B292-G010 
B337-G068 
B350-G162 


6.4±4.1 
9.2±3.3 
4.9±2.9 
9.3±2.3 


2.7±1.2 
2.6+1.9 
9.8±2.5 


5.9±3.4 
6.6±3.2 
12.4±4.7 


2.02+0.10 
1.00+0.10 
2.03+0.10 
1.99+0.10 


4.000+1.325 
0.003+0.000 
1.900+0.360 
1.278+0.148 



(") Ages from Puzia et al. (2005; P05); Beasley et al. (2005; Be05); Wang et al. (2010, W10); Fan et al. < f20T0l F10). B easley et al. ( [20051 derived 
two set of age estimates by comparing the same spectra with two different sets of SPSS models by Bruzual & Chariot (2003 , BC03) and Thomas 
et al. (12003 [ T03); we report both estimates in this table. Caldwell et al. ( 1201 It do not provide an age estimate for these clusters, but assign a 
conventional value of 14 Gyr to all of them. 



Table 3. Spectroscopic metallicities and [a/Fe] for the main targets . 



ID 


[Fe/H]" 


[ff/Fe] 




P05* 


G09 


Cll 


BOO 


F10 


P05 


B058-G119 r 
B292-G010 
B337-G068 
B350-G162 


-0.25+0.17 
-1.63+0.49 
-0.59+0.11 
-1.81+0.37 


-1.02+0.21 
-1.54+0.37 
-1.08+0.15 
-1.54+0.31 


-1.05+0.11 

-1.17+0.14 
-1.43+0.13 


-1. 45+0.24* 
-1.42+0.16 
-1.09+0.32 
-1.47+0.17 


-0.916+0.099 
-2.249+0.043 
-0.433+0.113 
-1.029+0.127 


-0.34+0.15 
0.11+0.31 
0.00+0.11 
0.38+0.28 



O [Fe/H] from Puzia etal. J20051 P05); Galled et al. ( 2009; G09); Caldwell et al. i20TT1 Cll); Barmby et al. ( 120001 BOO); Fan et al. J20T01 F10). 
(*) from Huch ra et a l. {1991). 

O Rich et al. ( 120051 provides a CMD-based estimate of [Fe/H] = -1.40 ± 0. 20 for the cluster B058. 
(*) Trasformed from [Z/H] with the formula of Thomas, Maraston & Bender ( 12003) . 



rected for the ACS geometrical distortions and the F435W and 
F606W were then cross-correlated to match common source^] 
This non-standard procedure prevented the DOLPHOT routine 
from calculating the required color term to perform the trans- 
formations from instrumental magnitudes to standard B,V mag- 
nitudes automatically; we made this final passage a-posteriori 
with the transformations by Sirianni et al. (2005). Below we will 
always consider only the photometry in the Johnson B,V pass- 
bands. Similarly to P09b, we retained in our final catalogs all the 
sources with valid magnitude measurements in both passbands, 
global quality flag = 1 (i.e., best measured stars), crowding pa- 
rameter < 0.5, x 2 < 2.5, and sharpness parameter between - 
0.3 and 0.3 (see Dolphin 2000b for details on the meaning of 
the various quality parameters). This selection cleans the sam- 
ple from the majority of spurious and/or bad measured sources 
without significant loss of information and was found to be ap- 
propriate for the whole data set. 

The limiting magnitudes of our photometry range from 
V~27.0 for the field of B058, which is projected onto a much 
more crowded background compared to the other clusters (see 
Figs.[T]and[2]i, to V~29.0 for the other fields. The internal photo- 
metric errors of individual measures are in the range 0.005 - 0.08 
mag for stars brighter than V=27 (see the error bars in Figs. [4] 
|5j and[6]). 

6 The match was performed with CataXcorr, a code aimed at cross- 
correlating catalogs and finding astrometric solutions, developed by P. 
Montegriffo at INAF - Osservatorio Astronomico di Bologna, and suc- 
cessfully used by our group for the past 10 years. 



2.1. Artificial stars experiments 

We performed artificial stars experiments (generating ~ 10 s fake 
stars per field) to study in detail the completeness of the samples 
as a function of magnitude, color, and distance from the clus- 
ter's center (i.e., as a function of crowding), using the automated 
procedure included in DOLPHOT (see P09b for details). 

The fields of B292, B336, and B350 all have completeness 
^80% for VS527, more than sufficient to detect a blue HB where 
it is present. On the other hand, the TO stars are located in 
a CMD region in which the completeness is less then 50%, 
which prevents us from reaching any suitable conclusion from 
this CMD feature. The field of B058 is projected onto the so- 
called 10 kpc ring, a site of ongoing star formation in the thin 
disk of M3 1 where the stellar density (and then the crowding) 
is much higher, conseguently the completeness is significantly 
worse (<;80% for V^24.5) and barely reaches the level of the 
red HB. 

2.2. Classification of secondary targets 

The CMD of the secondary target B049 obtained from these data 
is very similar to that already published and discussed in P09b, 
which was obtained from a different dataset; we refer the inter- 
ested reader to this previous publication for details. 

SK087B is the faintest object in our sample and is located 
in the dense disk field of B058. The CMD of the stars located 
within a radius of ~ 2.5" appears to be indistinguishable from 
that of the surrounding field and the same stars do not show 



4 



Perina et aL: HST/ACS CMDs of candidate intermediate-age M31 globular clusters. 



Table 4. Fundamental parameters and classification for addi- 
tional clusters and candidate clusters included in the images. 



ID 


RA (J2000) 


Dec (J2000) 


v„„ 


E(B-V)* 


type" 


22 
24 


B049 


00 41 45.6 


+40 49 53.7 


17.56 


0.16±0.02 


KD 


> 

26 


SK087B 


00 41 48.2 


+40 47 31.5 


18.83 




7(2) 




SK088B 


00 41 51.5 


+40 47 06.6 


16.98 




6(2) 


28 


SK089B 


00 41 54.7 


+40 47 14.6 


18.57 




2(2) 




B531 


00 36 14.7 


+40 57 56.6 


19.63* 




1 




B336 


00 40 47.6 


+42 08 41.4 


17.81 


0.42±0.04 


KD 


22 



V,-„ , from the RBC. * from Fan et al.( 2008 2010). Caldwell et 
al. f20TTt reported E(B-V)= 0.32 and 0.05 for B049 and B336, 
respectively. ": RBC Classification flag: 1- confirmed cluster; 2- gc 
candidate; 6- star; 7- asterism. The classification reported in the RBC 
before this study are enclosed in parentheses. * See Sect. |3.2.6| 
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any over-density with respect to the field. For these reasons re- 
classified this object as an asterism. The high-resolution ACS 
images clearly reveals that SK088B is not a genuine cluster but a 
bright star, instead. SK089B is imaged only on the frame F606W, 
so we cannot provide its CMD, however, it appears to be an OB 
association in the proximity of its dust cocoon. We maintained 
the classification of this object. 



3. The CMDs 

Below we present and discuss the CMDs of the four main tar- 
get clusters, of the secondary target cluster B336 and of the 
newly detected cluster B531. The individual CMDs are shown 
in Figs. |4j [5] (target clusters), and [6] (secondary clusters). The 
cluster CMDs shown in the left panels sample the stellar popula- 
tion within an annulus around the cluster center. The inner limit 
of the annulus is set by the crowding level at the center, which 
prevents us from performing useful photometry in the most cen- 
tral region of the cluster, the outer limit is set by the need to 
avoid strong contamination by the surrounding field population. 
The inner and outer radii of the adopted annuli are indicated for 
each cluster. In the right panels we plot the CMD of the field 
adjacent to each cluster. The field population is sampled over a 
wide region outside the outer limit of the cluster and the area of 
the considered field is reported in units of the area of the annulus 
enclosing the cluster stars. 

In the left panels of Figs. [4] a nd [5] w e over-plotted the 
isochrones from the Girardi et al.( [2002 ) set on the cluster 
CMDs, taking the age and the metallicity estimated by P05 for 
these clusters. The isochrones have been shifted to the same dis- 
tance and reddening of the clusters adopting the values derived in 
Sect. |3 - 1 1 and listed in Table [6] Clearly, the CMDs are not deep 
enough to reliably test the region of the upper MS and MSTO 
that is expected for populations of that age. The handful of stars 
lying around the expected sub-giant branch (SGB) level are fully 
consistent with being largely caused by field contamination (see 
Fig. [8J. In any case, we can obtain strong and meaningful lower 
limits to the age of three of the four main targets (excluding 
B58) by searching for the oldest isochrone whose SGB clearly 
is above the observed candidate SGB stars (see Fig. 8] below). 
The resulting limits, (age > 6.3 Gyr for B292, age > 4.5 Gyr for 



Fig. 4. Left panels: CMD of the stars enclosed in the reported 
annulus for the target clusters B292 and B350. Overplotted on 
the CMD there are the 80% and 50% completeness levels (dotted 
lines) and an isochrone from the Girardi et al.( 2002) dataset with 
age and metallicity corresponding to the spectroscopic estimates 
by P05 (see Tableland [3}; the red giant branch is plotted with 
a thick continuous line and the asymptotic giant branch with a 
thin dashed line. The reddening and distance moduli are taken 
from Fig. [7] below. Right panels: CMD of the field stars in the 
neighborhood of the cluster. The area of the considered field is 
reported in units of the area of the annulus shown in the left 
panels (A). 




10 12 
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Fig. 5. Same as [4] for the other two target clusters B058 and 
B337. 
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Fig. 6. Same as Fig.^but for the clusters B336 and B531. The 
CMD of B049 is very similar to that obtained by P09b and is not 
shown here. 



B337, and age > 5.6 Gyr for B350) are sufficient to rule out the 
estimates from integrated photometry by W10 and F10, which 
point to ages significantly younger than P05 (see Table |2j. On 
the other hand, we have to rely on evolved stars to constrain the 
age of our clusters in the range invoked by P05. 

In most CMDs the evolved population of the cluster can 
be clearly distinguished from the field (see also the de- 
contaminated CMD presented in Sect. |3.1| l. In particular, B292 
and B350 display an obvious blue HB that cannot be accounted 
for by field contamination (see Fig. [7] and Table |7J. As said, 
this is a characteristic feature of very old metal-poor GCs; its 
mere presence strongly argues against ages younger than 10-11 
Gyr. For the same clusters, the over-plotted isochrones predict 
the (possible) presence of asymptotic giant branch stars brighter 
than the RGB tip (dashed line) while no such stars are observed, 
providing further support to the idea that these clusters are in- 
deed older than the considered models. It is worth to stress here 
that for B292 and B3 50 an a ge of 12.5 Gyr (i.e., typical of old 
MW GCs, Dotter et al. 2010) is within Icr from the age estimates 
by P05, hence fully compatible with them (see Table |2j. It may 
be interesting to investigate why the accurate Lick indices analy- 
sis performed by those authors singled out these clusters as being 
younger than the bulk of the others. We will attempt to achieve 
some insight into this question in Sect. [5] On the other hand, 
even an age ^10 Gyr would strongly disagree with the values 
provided by W10 and F10, suggesting that these estimates are 
not only systematically too low, but appear to have unrealisti- 
cally low associated errors (typically lower than 1 Gyr; Table [2j. 

3. 1 . Determination of the cluster parameters 

Before proceeding with the analysis of the cluster properties, 
we applied the field decontamination procedure described in 
Bellazzini et al. dl999l l. This method is based on a clipping rou- 
tine which, making use of the local density on the CMDs of the 



Table 5. Reference grid of template Galactic globular clusters. 



ID [Fe/H] E(B - V) (m - M) 







mag 


mag 


NGC7078 (M15) 


-2.16 


0.10 


15.37 


NGC5824 


-1.87 


0.13 


17.93 


NGC6205 (M13) 


-1.65 


0.02 


14.48 


NGC5904 (M5) 


-1.40 


0.03 


14.46 


NGC6723 


-1.12 


0.05 


14.85 


47Tuc 


-0.71 


0.04 


13.37 


NGC6624 


-0.35 


0.28 


15.36 



NOTES: Metallicities are from Zinn (1985); all other parameters are 
from Harris (1996) (online update 2003). B,V photometry is from 
Piotto et al. (2002). 




0.5 1 1.5 2 0.5 1 1.5 2 



(B-V) (B-V) 

Fig. 7. CMDs of the clusters after statistical decontamination 
from field stars with a superposed grid of HB and RGB tem- 
plates. The ridge lines that provide the best fit to the observed 
HB and bracket most of the RGB population are plotted as heav- 
ier lines to highlight them. The adopted values for the metallic- 
ity, reddening, and distance modulus are reported in the upper 
left panel of each plot. 

field and of the cluster, computes the probability that a given 
star is a member of the cluster and retains or rejects stars from 
the cluster CMD according to that (see P09b for further details). 

Assuming that the clusters are classical old globulars, we 
derived a simultaneous estimate of the distance modulus, red- 
dening, and metallicity by comparing them with a set of CMD 
templates of well studied Galactic GCs, exactly in the same way 
as in P09b. We searched for the set of parameters ((m - M)o , 
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Table 6. Cluster parameters from the analysis of the CMD. 



ID 


E(B-V) 


(m - M) 


[Fe/H] 




mag 


mag 




B292-G010 


0.12 


24.43 


-1.80±0.25 


B337-G068 


0.06 


24.43 


-1.30±0.25 


B350-G162 


0.08 


24.45 


-1.50±0.25 


B336-G067 


0.10 


24.45 


-1.8 ±0.25 


B531 


0.10 


24.43 


>-0.42 



Table 7. Star counts in CMD boxes enclosing the BHB for the 
three clusters B292, B350, and B336. 



ID 


N GC 
obs 


fffl-U 
obs 


N GC 
1 Y exp 


cr 


AreaQc 


B292-G010 


30 


25 


1.05 ±0.21 


5.3 


0.042 


B350-G162 


37 


76 


2.36 ± 0.27 


5.7 


0.031 


B336-G067 


23 


21 


0.55 ±0.12 


4.7 


0.026 



is the number of cluster stars expected in a given box from the 
field population, computed by rescaling the observed counts in the 
sampled field by the ratio between the sampled GC and field areas (last 
column). In the fifth column we list the background-subtracted star 
counts in the GCs sampled areas in unit of a. 



E(B-V) and [Fe/H]) producing the best match between the ob- 
served RGBs and HBs and the ridge lines of the template clus- 
ters in the absolute plane, given the direction of the reddening 
ve ctor A B = 4.145£(fi - V), A v = 3AE(B - V) (Schlegel et 
al. 119981 ). The best match was judged by eye guided by (exten- 
sive) experience. Color and magnitude shifts are applied itera- 
tively until a satisfactory match with any RGB and HB template 
is found. From these shifts we obtained estimates of the redden- 
ing and distance, while the metallicity was estimated by inter- 
polation between the two RGB ridge lines bracketing the ob- 
served RGB locus. As starting values for the iterative procedure 
we used E(B-V) = 0.08 for the foreground reddening (Barmby 
et al. 2007; Burstein & Heiles 1984) and the distance modu- 
lus (m - M)o = 24.47 (McConnachie et al. |20D31 >. The ridge 
lines of the reference GGCs were derived from publicly avail- 
able photometry (Piotto et al. 2002) and were shifted to the ab- 
solute reference frame by correcting for reddening and distance 
using the values listed in Table [5] These reference GGCs were 
chosen to provide a sufficiently fine and regular sampling over a 
wide enough range of metallicities for a correct characterization 
of the target GCs. 

In Fig. [7] we show the field-decontaminated CMDs with the 
reference grid of GGC ridge lines over-plotted. The values of 
metallicity, reddening, and distance corresponding to the best 
match are also reported in each individual panel, as well as 
in Table [6j the typical uncertainty on the distance modulus is 
0.2 mag, 0.04 mag in E(B-V), and 0.25 dex in metallicity. In 
Sect. |3.2| we briefly discuss each individual cluster. 



3.2. Results for individual clusters 

3.2.1. B292=G010 

The cluster B292 lies at a fairly large projected distance from the 
center of M31, in a low -density region where the contamination 
by the field stars is very low. Its CMD (Fig.|4j top-left panel), is 
characterized by a steep RGB and a well defined and populated 
BHB that is obviously associated with the cluster: star counts 
within a box enclosing the feature indicate an overabundance of 
BHB stars in the cluster with respect to the surrounding field (see 
Table [7j at 5 cr level. 

The best match with the Galactic RGB and HB ridge lines 
is obtained by assuming a value of reddening E(B - V) = 0.12 
mag and a distance modulus of (m - M)q = 24.43. With these 
assumptions, the cluster's RGB fall between the ridge lines of 
NGC5824 and M13, indicating that the metallicity of B292 is 
[Fe/H] = -1.8 ± 0.25. This value agrees well with the spec- 
troscopic metallicity estimates by P05 and Galleti et al. (2009 
hereafter G09), while it is just marginally more metal-poor than 
what was found by Barmby et al. (2000 , hereafter BOO), and 0.45 
dex more metal-rich than the value obtained by F10 from inte- 
grated photometry. The mismatch in color between the vertical 
parts of observed and template BHB may suggest a solution with 
a higher reddening and, consequently, a lower metallicity; this 
would not change the main conclusions of the present analy- 
sis significantly. We preferred the solution presented in Fig. [7] 
which gives greater weight to the five stars at intermediate color 
(0.1 < (B - V) < 0.6, around M v =s 0.6), because, if inter- 
preted as genuine cluster members, they provide a constraint on 
the vertical direction (distance) that is otherwise lacking. The 
best match of the color of the vertical portion of the HB can be 
obtained by adopting E(B-V)=0.19, and (m - M) = 24.3. This 
solution appears unlikely because it would imply that this clus- 
ter is located several kpc above (and in front of) the M3 1 disk 
while suffering from a degree of extinction quite typical of clus- 
ters embedded in the disk (see Perina et al. 2010 ). A small-scale 
Galactic cloud on this line of sight seems the only possibility to 
make this solution viable. In any case, this solution also yields a 
low metallicity ([Fe/H] =s —1.9) and the old age implied by the 
presence of BHB. 

3.2.2. B350=G162 

The CMD of B350 (Fig.|4| bottom-left panel) shows a well pop- 
ulated blue HB clearly associated with the cluster (see Table [7J 
and Fig. [8]). The RGB is steep, which indicates a low metal con- 
tent. 

The best match of the RGB and HB features with the cor- 
responding Galactic ridge lines is obtained by assuming a dis- 
tance modulus of (m - M)o = 24.45 and maintaining the adopted 
starting value of E(B - V) = 0.08 mag for the reddening, in 
good agreement with what was reported by F10 (Table [l}. With 
these assumptions, the cluster RGB fall between the ridge lines 
of M13 and M5, indicating for B350 a metallicity of [Fe/H] = 
-1.5 + 0.25. This is compatible with most of the values avail- 
able in literature, except for the marginal discrepancy with the 
estimate by F10 from integrated photometry (see Table[3]). 

3.2.3. B058=G119 

Owing to the high degree of crowding, the CMD of B058 is not 
of sufficient quality to obtain a truly reliable solution. The CMD 
previously obtained by Rich et al.( 2005 ) from WFPC2 data 
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Fig. 8. CMDs of the candidate intermediate-age clusters (stars 
selected in the same annuli as in Figs. [4] and [5]) before (left 
panels) and after (right panels) the statistical decontamination 
from field stars. We superposed on the CMDs the oldest-age 
isochrones whose upper MS and SGB trace the blue/bright edge 
of the observed stars to set a strong lower limit to the age of the 
clusters (metallicity, reddening, and distances are the same as 
those adopted in Figs.[4]and[5]l. The age of the adopted isochrone 
is reported as a label in the upper-right corner of the left pan- 
els. Note, however, that there are no such faint stars in the de- 
contaminated CMD, suggesting that they are not associated to 
the clusters. 



in the filters F555W (V) and F814W (I) with longer exposure 
times (5300 s and 5400 s, respectively) seems easier to inter- 
pret, because it appears to be slightly deeper and the RGB shape 
is straighter in the I-band. The solution presented here must be 
considered as only tentative, we suggest to keep the parameters 
by Rich et al. (120051) as the best set. 

The 6.4 Gyr isochrone shown in Fig. [5] seems fully com- 
patible with the observed CMD, but an old isochrone would 
fit equally well. The estimate by P05 had a large uncertainty 
(6.4 + 4.1 Gyr) and the CMD test is completely inconclusive. 
Ages as low as 4 Gyr or 2 Gyr, as proposed by F10 and W10, 
respectively, are compatible with the observed CMD as well. 

3.2.4. B337=G068 

B337 has a very clean CMD characterized by a steep RGB and 
a well populated red HB (see Fig. [5j. The decontamination has 
little effect on the CMD of this cluster, which is positioned in a 
quite low-density region (see Figs. [7] and [8}. 

The best match of these features with the corresponding 
ridge lines is obtained by assuming a value of reddening of E(B— 



V) = 0.06 mag (in excellent agreement with F10) and a distance 
modulus of (m - M)q = 24.43. With these assumptions, the clus- 
ter RGB lies between the ridge lines of M5 and NGC6723, in- 
dicating a metallicity of B337 of [Fe/H] = -1.3 ± 0.25. This 
value agree with the estimates by G09 and BOO, while it is sig- 
nificantly lower than what was obtained by P05 and F10. (see 
Table |3j. The CMD appears to be incompatible with this high 
metallicity, if genuinely old. 

It is impossible in this case as well to draw a firm conclu- 
sion on the age of the cluster. The red HB is bluer than that of 
the field, compatible with old metal rich population a la 47 Tuc. 
Three stars lying very close to the BHB templates are preserved 
in the de-contaminated CMD, but they are too few to be un- 
equivocally attributed to the cluster. As said above, the faint end 
of the CMD seems already incompatible with the formal best-fit 
age proposed by P05 (4.9 Gyr; see Fig. |8j, but the tension re- 
laxes if the uncertainties are taken into account (the same is true 
for the estimate of Beasley et al. 2005 ). On the other hand, the 
values proposed by F10 (1.9 + 0.4 Gyr) and W10 (2.0 ± 0.1 Gyr) 
are clearly too low to be compatible with the CMD. 

3.2.5. B336=G067 

B336 is a secondary target falling in the field of B337. The 
CMD of B336 presented here shows a populous intermediate- 
blue HB population, implying an old age for the cluster (see 
Table fl). This is also the only cluster in our sample that dis- 
plays a population of candidate RR Lyrae variables in signifi- 
cant excess with respect to the surrounding field, another indi- 
cation of old age (Sect. [4j. The best fit to the decontaminated 
CMD with the Galactic templates is obtained assuming a dis- 
tance modulus of 24.45 and a E(B - V) = 0.10. The RGB of 
B336 lies between the ridge lines of NGC5824 and M13, in- 
dicating a metallicity [Fe/H] = -1.8 ± 0.25 (see Fig. |7J. Fan 
et al. (120101) estimates from integrated photometry an age of 
1 .015 ±0. 167 Gyr, a value incompatible with the observed CMD. 
The authors obtained a metallicity of [Fe/H] - 2.221 ± 0.085 
in marginal disagreement with the value estimated here, and a 
reddening E(B - V) = 0.420 ± 0.040. The adoption of this high 
reddening value would shift all RGB and most BHB cluster stars 
to the blue of the bluest template, leading to an obvious incon- 
sistency. Caldwell et al. d20TTT > report [Fe/H] - 2.49 ± 0.59 and 
E(B-V) = 0.05, confirming that the cluster is not very reddened 
and is remarkably metal-poor. 



3.2.6. B531 

B531 is a newly detected cluster, found in the field of the target 
B292. Its CMD in Fig. [3] shows a red HB and an RGB quite 
tilted and redder then the reddest available Galactic template 
(NGC6624). The fitting procedure of the decontaminated CMD 
with the Galactic templates provided the best result adopting a 
distance modulus of 24.43 and E(B - V) = 0.10 mag (similar 
to what was found for B292 in the same field), indicating B531 
as a very metal-rich cluster with [Fe/H] > -0.42 (see Fig.|7i. 
A solution with a higher reddening would lead to a lower meta - 
licity estimate, but the match of the red HB level would be lost. 
Moreover, this line of sight is far from the visible disk of M31 
and is unlikely to be associated with high values of intrinsic ex- 
tinction. B531 lies at a projected distance R p 17 kpc from 
the center of M31; the projection along the galaxy minor axis 
is Y ~ 11 kpc. In P09b it has been noted that metal-rich clus- 
ters located so far away from the galaxy center (and "above" 
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Fig. 9. Process of detection and selection of candidate RR Lyrae 
is illustrated using B336 as an example. Upper panel: the stan- 
dard deviation of the five individual measures of F606W is 
plotted as a function of the average F606W. The candidates 
RR Lyrae were selected among stars with 24.9 < F6Q6W < 
25.4, as stars with rms deviations larger then 3 cr (open stars). 
The continuous line indicates the 3 <x level and the dotted lines 
the 2 and 1 <x levels. Lower panels: Light curves of the candi- 
dates RR Lyrae detected in the upper panel and labeled with 1, 
2, and 3, taken as examples. Candidates 1 and 2 show coher- 
ent patterns of variations at all sampled epochs, compatible with 
genuine RR Lyrae variables. Candidate 3 is a classical case of 
spurious detection: F606W is constant at all epochs except one, 
likely corresponding to the hit of a cosmic ray over the PSF area 
of the star in the individual image at that epoch. 



the galaxy disk) are quite rare and may trace substructures as- 
sociated to past accretion events in the M3 1 halo (see the case 
of B407, discussed in P09b, and also Mackey et al. l20T0t . If the 
high-metallicity of B53 1 will be confirmed, this cluster would be 
an interesting candidate "accreted" GC. A spectroscopic study, 
aimed at estimating the metallicity and the velocity of this clus- 
ter would be very valuable to follow-up this hypothesis. The in- 
tegrated color of B531 (B - V)o = 0.81 is very similar to that of 
47 Tuc [(B - V) = 0.84, from Harris [19961 . thus supporting a 
(relatively) high metal content. 

Integrated B,V, I magnitudes were computed for this clus- 
ter by transforming aperture magnitudes obtained from Sloan 
Digital Sky Survey (Abazajan et al. 120091) calibrated images, as 
described in detail in Perina et al. 2009a for VdBO. Adopting an 
aperture of radius 11 px^ 4.3", we obtained g = 20.08 + 0.04, 
r = 19.32 ± 0.04, and i = 18.99 ± 0.04. Lupton (2005) transfor- 
mations^] were used to convert from g,r,i to B,V, I magnitudes. 



www. sdss .org/dr4/algorithms/sdssUB VRITransform.html 



4. Detection of RR Lyrae variables 

Another possibility to classify a cluster as very old is to identify a 
population of RR Lyrae belonging to it. The considered observa- 
tions were not planned for a search of variables and are far from 
optimal for this purpose. For each observed field, the observa- 
tions in the two filters were acquired at a distance of two months, 
therefore they are not expected to sample the same phase range 
of any light curve. The four consecutive F435W images (per 
field) span an interval of ~ 3.2 hours (i.e., ~ \ of the typical 
period of a RR Lyrae); the six consecutive F606W images (per 
field) span an interval of ~ 4.0 hours (i.e., ~ i of the typical 
period of a RR Lyrae). Therefore, even if the amplitude of the 
variation is expected to be larger in F435W, the F606W observa- 
tions are likely more efficient in detecting RR Lyrae because of 
the larger number of exposures and larger phase coverage. 

With this observational material we can hope to detect ab- 
type RR Lyrae variables with periods longer then ~ 5 h only 
if they are caught during the rapid rise to maximum light or in 
the descending branch of their light curve, while they would go 
unnoticed if they are observed near their maximum/minimum. 
Moreover, because the exposures in the two filters are not con- 
secutive, the detection of variability in both passbands is only 
a serendipitous event. Similar limitations have been encoun- 
tered also by Clementini et al. (2001), who for the first time, 
searched for RR Lyrae variables in a sample of M31 GCs us- 
ing archive HST Wide Field Planetary Camera 2 (WFPC2) data. 
That study, however, proved the possibility of a successful de- 
tection of likely candidates RR Lyrae also with data not opti- 
mized for this type of investigation. It is important to stress that 
the present analysis is not intended to obtain a complete census 
of the RR Lyrae population in the considered clusters but, in- 
stead, is aimed at verifying whether a significant population of 
RR Lyrae can be identified in some of them. 

To identify RR Lyrae candidates, we adopted the following 
criteria, which are ilustrated in Fig. [9]: 

1. We searched for variability only in stars lying within + \mag 
of the HB level as estimated in Fig. [7] 

2. We considered objects with rms deviations in their estimate 
of the magnitude larger than or equal to three times the av- 
erage deviation at their luminosity level. Fig.[9]shows an ex- 
ample of application to the cluster B336: the selection pro- 
cedure identifies 1 1 candidate variables (open star symbols) 
that satisfy this criterion. 

3. For the selected stars the sequences of individual measure- 
ments were then displayed as a function of the epoch, search- 
ing for trends in the light curves that were compatible with 
RR Lyrae type variations. Fig. 10 shows as an example 
the light curves of the three candidate variables labeled in 
Fig. [9] The stars labeled 1 and 2 show a coherent pattern 
of variations in all sampled epochs, compatible with gen- 
uine RR Lyrae variables. On the other hand, candidate 3 is a 
classical case of spurious detection: F606W is constant at all 
epochs except one, likely corresponding to the hit of a cos- 
mic ray over the PSF area of the star in the individual image 
at that epoch. 

4. We identified a star as a candidate RR Lyrae variable if it 
satisfied all three criteria, at least in one of the two passbands. 

As a result of this analysis, we found 1 1 (B336), 3 (B058), 2 
(B292, B350), 1 (B531) and (B337) candidate RR Lyrae stars 
in our six clusters, respectively. Their positions in the CMDs are 
shown in Fig. 10 We searched for candidate RR Lyrae variables 
also in comparison fields surrounding the clusters to quantify the 
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Fig. 11. Examples of M31 globular clusters classified as Blue 
Horizontal Branch (BHB, upper panels) or Red Horizontal 
Branch (RHB, lower panels). All the CMDs are from Rich et 
al. ( |2005T l except B407, which is from P09b. 



Fig. 10. Candidate RR Lyrae are overplotted as open stars on the 
cluster CMD. 



probability that the variables detected in the cluster can be ex- 
plained as merely due to the field population. Given the observed 
frequency of candidate variables in the surrounding fields, the 
number of field RR Lyrae expected to fall on the cluster area 
are 0, for the clusters B336, B292, B531, B337, 1, for B350, 
and 2 for B058. Hence, the only cluster for which a significant 
excess of RR Lyrae is found with respect to the expectations is 
B336, providing further evidence supporting a very old age for 
this system. The very red colors of most of the candidates in 
B058, B350, and B531 also cast some doubt on their classifica- 
tion as RR Lyrae. Note that the lack of detection of a significant 
excess of candidate RR Lyrae in the other clusters cannot be 
taken as evidence that no such variable is present there: our data 
are clearly not sufficient to reach this conclusion. 



5. Summary and discussion 

The main aim of this study was to verify whether the interme- 
diate age estimate obtained from integrated spectroscopy and/or 
photometry by several authors (Table |2|l for the M3 1 clusters 
B058, B292, B337, and B350 was supported by the CMD ob- 
tained from the new HST-ACS observations. The final answer is 
completely inconclusive for B058 and uncertain for B337, even 
if in the last case there is marginal evidence for an older age com- 
pared to the dating from integrated light. On the other hand, both 
B292 and B350 show a significant and unequivocal population 
of blue horizontal branch stars, indicating that they are as old as 
the classical Galactic GCs (ages older than 1 1 Gyr; see Dotter et 
al. 2010, for a recent analysis and discussion). It should be noted 
that the age estimates provided by P05 are consistent within 1 cr 
uncertainties with this conclusion. However, these clusters were 



singled out as likely younger than the bulk of M3 1 clusters at 
similar metallicity and this classification is clearly not confirmed 
here. On the other hand, ages as young as claimed by W10 and 
F10 (all < 4 Gyr) are clearly ruled out by the observed CMDs 
for all main targets of the present study except (possibly) B058. 
The comparison with the results by Be05 significantly depends 
on the considered set of age estimates provided by these authors 
(see Table[2]l. It may be interesting to note that except for B350, 
their estimates based on the BC03 SPSS models are not com- 
patible (too young) with the observed CMDs, even taking into 
account the reported uncertainties; on the other hand, the esti- 
mates obtained from the T03 models are similar to those by P05, 
and, in particular, the old age reported for B350 fully agrees with 
our results. 

The sensitivity of colors and/or spectral indices generally 
used as age diagnostics (like, for example, H/3) to the morphol- 
ogy of the horizontal branch is known for long time: for example, 
it was already taken into account in the models of the integrated 
light of SSPs by Buzzoni dl989l l. Bright and hot BHB stars typi- 
cal of an old population may mimic the effect of a hotter MSTO, 
i.e., of a younger age (see Lee et al. l2000l for a thorough discus- 
sion, and Percival & Salaris l2010l for a very recent investigation 
in the context of SPSS modeling). It is tempting to identify this 
as the main reason at the origin of the erroneous classification of 
B292 and B350 as intermediate-age clusters (Rey et al. 120071 . 

To gain a deeper insight into this hypothesis, we considered 
all M31 GCs for which a safe classification of their HB mor- 
phology can be obtained from their CMDs, taken from different 
sources (see Sect. [T] for an exhaustive list). We simply divided 
the clusters into BHB (blue squares) and RHB (red circles) from 
the inspection of their CMD; our classification criterion is illus- 
trated in Fig. 1 1 We ended up with a sample of 33 clusters that 
also have estimates of the H/3 and Mg2 indices reported in the 
RBC (see G09). The classical age-diagnostic plot Mg2 vs. H/3 
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Perina et aL: HST/ACS CMDs of candidate intermediate-age M31 globular clusters. 



X 



X 




Fig. 12. Galactic (lower panel) and M3 1 (upper panels) globu- 
lar clusters whose HB morphology can be classified from their 
CMD in the Mg2 vs. H/3 plane. The average 1 cr error-bars are 
plotted in the upper part of the panels. The over-plotted grid of 
models is from the recent set by Thomas et al. (2010): contin- 
uous lines are iso-age loci of age 3, 5, 8, 12 and 15 Gyr as la- 
beled in the lower panel; dotted lines are iso-metallicity loci for 
[Z/H]=-2.25, -1.35, -0.33, 0.0 and +0.35, as labeled on top of the 
lower panel. The adopted models are for SSPs with [ar/Fe]=+0.3. 
Blue squares are clusters classified as HB-blue, red-circles are 
clusters classified as HB-red. The main targets of this study are 
represented with open symbols and labeled. We also labeled 
three M31 clusters lying above the age=5 Gyr line (B010, B293 
and B008) and two Galactic clusters lying on top or above the 
age=8 Gyr line (NGC5024 and NGC3201). The arrows on the 
points for B337 and B008 show the position implied by the new 
estimates of H/3 (Galleti, private communication). 



for these clusters is shown in the upper panel of Fig. 12 a grid of 
SSP models from the most recent set by Thomas et al. (1201 Oi l is 
over-plotted as a reference. None of the conclusions drawn from 
this plot is significantly dependent on the actual set of models 
adopted (see, for example, Fig. [13] below). In the lower panel 
of Fig. [T2]the same diagram is shown for the Galactic GCs. The 
H/3 and Mg2 indices are taken from G09. In this case we classi- 
fied as BHB clusters with B 1 R * V > 0.0 and as RHB those with 
B 1 R 1 V < 0.0, where the values of the classical HB morphology 
parametei *] gf^fy are taken from Harris d 19961 1. 

In addition to the obvious correlation between the HB mor- 
phology and the metallicity of the clusters^ there are several 
other interesting features in Fig. 
mented on in the present context. 

8 Where B, R, V are the number of stars to the blue, to the red, and 
within the RR Lyrae instability strip, see Lee, Demarque & Zinn ( 1990 1. 

9 Metallicity is the first parameter driving the HB morphology, see 
Sa ndage & Wildey ( 11967b . Fusi Pecci & Bellazzini d!997t . Dotter et 
al. J2010t . and reference therein. 



12 that deserve to be com- 



1. According to Fig. [12] B292 and B350 appear younger than 
8 Gyr, in agreement with P05. These are not the most ex- 
treme cases, because B0 10 and B293 have even larger H/3 
values, implying ages as young as ~ 3 Gyr and ~ 2 Gyr, re- 
spectively, if the comparisons with models shown in Fig. [12] 
are taken at face value. Nevertheless, also these clusters, as 
B292 and B350 (and also B311, discussed in Sect. [I), dis- 
play an obvious BHB population (see, e.g., Fig. [TT] for the 
case of B293) indicating that they are, in fact, much older 
than this. 

2. In this context it is interesting to note that the well stud- 
ied BHB Galactic clusters NGC5024 and NGC3201 fall, in 
the same diagram, on top and just above the age=8 Gyr lo- 
cus, respectively. We have reliable age estimates for these 
clusters based on exquisite photometry that reaches well 
below the MSTO. In the most recent thorough analysis, 
based on an homogeneous set of HST data, Dotter et 
al. d20T0l found an age=12.0+0.75 Gyr for NGC3201, and 
an age=13.25+0.50 Gyr for NGC5024, in an age scale in 
which the oldest GC has age= 13.50+0.50 Gyr. This shows 
quite conclusively that even in the best possible case, an age 
classification based on the classical spectral diagnostic H/3 is 
not at all able to unambiguously distinguish between old and 
intermediate-age clusters, because 13 Gyr old BHB clusters 
can easily have H/3 ~ 3A, i.e., the same as ~ 5 - 8 Gyr old 
clusters (see Lee et al. 2000). In this sense, it appears reason- 
able the limit H/3 > 3. 5 A adopted by Fusi Pecci et al. (120051 
and in the RBC) to identify genuinely young clusters, even 
if this a conservative choice may still not be completely safe 
from false positives, like the case of B008, described below. 

3. Returning to M31 BHB clusters, the situation shown in 
Fig. [12] may appear even worse. However, the larger spread 
of the distribution of M31 points with respect to MW ones is 
suggestive of the effect of larger observational errors, which 
would not be unexpected given the large difference in dis- 
tance between the two sets of GCs. Even if average formal 
errors provided in the literature are similar in the two cases 
(see the error-bars in Fig. [L2| , the fact that the M31 BHB 
clusters show similar spreads above and below the 12-15 Gyr 
models is consistent with our hypothesis. 

4. Further support of this view was provided by new spectra 
that were recently obtained for two RHB clusters that ap- 
peared to have an anomalously high H/3, B337 (one of our 
main targets) and B008 (lying straight on the age=l Gyr lo- 
cus j^] The newly obtained estimates of H/3 for these clusters 
are much lower than those found in the literature (see the ar- 
rows in Fig.[L2| and fully agree with ages ^ 12 Gyr. This new 
finding suggest that large errors in Lick indices taken from 
the literature may not be such rare occurrences for M3 1 GCs 
as one might expect. 

To illustrate the effects of the HB morphology on Hp, we 
provide in Fig. 13 the same comparison as in Fig. 12 but using 
models from the set by Buzzoni (1989). The overall results are 
similar to those obtained in Fig. 12 hence we do not repeat the 
discussion. The key point displayed in Fig. 13 is that if the ob- 



served range of HB morphology is taken into account, the H/3 

10 S. Galleti, private communication. The spectra were obtained at the 
Cassini Telescope, in Loiano, using the same set-up and data-reduction 
procedure described in Galleti et al. ( 2005 , 2009 1 The H/? index was 
measured as in Galleti et al. d2009b . It is found H/3 = 1.67 ± 0.39 A 
instead of H/3 = 3.67 ± 0.30 A for B008, and H/3 = 2.37 ± 0.36 A 
instead of H/3 = 3.23 ± 0.07 A for B337. 
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Perina et aL: HST/ACS CMDs of candidate intermediate-age M31 globular clusters. 
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Fig. 13. Same as Fig. 
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but in this case the comparison is 
with models from the Buzzoni (1989) set, to illustrate the ef- 
fects of the HB morphology on H/3. The lower continuous line 
is an age=15 Gyr model that splits in two branches around 
Mg2 = 0.12 (i.e., in the low metallicity regime) according to 
the HB morphology associated to the model, i.e., red (RHB) or 
intermediate (IHB). The latter was chosen because it provides 
H/3 values (slightly) higher than the BHB case. The dashed line 
is an age=15 Gyr model with enhanced oxygen abundance to 
show the effect of a-enhancement. 

values of all Galactic GCs can be reconciled with the old ages 
estimated from the CMDs by Dotter et al. (1201 Oi l. 

In conclusion, a convincing case has still to be provided for 
the presence of an intermediate-age bright cluster in M31 (see 
Caldwell et al. 120111 for a closely similar conclusion, reached 
with a different observational approach). The experiment per- 
formed here suggests that HST observations of M31 clusters 
with "reasonable" exposure times (^2-3 HST orbits) can re- 
liably probe only ages lower than 4-5 Gyr, while they would 
be inconclusive for older ages (except in the case of detection of 
BHB stars). On the other hand, the results of the present analysis 
clearly indicate that while the analysis of classical Lick indices 
remain a valuable technique to identify candidate intermediate- 
age GCs in external galaxies (including M31), they do not pro- 
vide sufficient evidence to firmly establish the case of ages in 
the range ~ 5 - 9, but need independent confirmation from other 
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